Heat transport through a Josephson junction 
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We discuss heat transport through a Josephson tunnel junction under various bias conditions. We 
first derive the formula for the cooling power of the junction valid for arbitrary time dependence of 
the Josephson phase. Combining it with the classical equation of motion for the phase, we find the 
time averaged cooling power as a function of bias current or bias voltage. We also find the noise 
of the heat current and, more generally, the full counting statistics of the heat transport through 
the junction. We separately consider the metastable superconducting branch of the current-voltage 
characteristics allowing quantum fluctuations of the phase in this case. This regime is experimentally 
attractive since the junction has low power dissipation, low impedance and therefore may be used 
as a sensitive detector. 



I. INTRODUCTION 



Thermoelectric effects in metallic nanostructures are 
subject of intense experimental and theoretical research 1 . 
In this paper we consider a particular example of such a 
structure - a Josephson junction between two supercon- 
ductors. Although Josephson junctions are very impor- 
tant building blocks of many low temperature devices, 
their thermoelectric properties are still not fully under- 
stood and there exist only a few theoretical studies of this 
issue in the literature. Guttman et a£ have derived the 
heat current through a Josephson junction with different 
temperatures of the leads and biased below critical cur- 
rent so that the Josephson phase was fixed. Subsequently, 
Zhao et have corrected a sign error in their result 
and extended it to the case of arbitrary transparency of 
the barrier between the superconducting leads. Frank 
and Krech£ have considered a voltage biased Joseph- 
son junction and demonstrated that the superconduct- 
ing lead with a smaller gap gets cooled in a certain range 
of bias voltages. They have derived an expression for 
the cooling power of the junction, which was a straight- 
forward generalization of the well known result for the 
normal metal - insulator - superconductor (NIS) tunnel 
junction 6 ' 7 . To the best of our knowledge, the systematic 
experimental studies of the thermoelectric properties of 
Josephson junctions are still lacking. 

In this paper we develop a general approach to this 
problem, which covers the two regimes discussed above 
as limiting cases. First, we allow the Josephson phase 
to depend on time in an arbitrary way and derive the 
expression for the heat current, or cooling power, under 
these conditions. Next, we determine the actual time 
dependence of the phase by solving the differential equa- 
tion of motion, which describes phase dynamics under 
given bias conditions and in the presence of noise of the 
environment^.. More specifically, we use the model of 
rcsistively and capacitively shunted Josephson junction 
(RCSJ model) with the shunt resistance much smaller 
than both the resistance quantum and the resistance of 
the junction itself (see Fig. [T]). Afterwards, we find the 
cooling power averaged over time and noise which can be 



measured in an experiment. 

At bias currents below the critical one we use a more 
elaborate technique treating the Josephson phase as a 
quantum operator. It is necessary to ensure the detailed 
balance between forward and backward rates of the quasi- 
particle tunneling. In this limit we study the dependence 
of the cooling power on the bias current and on the differ- 
ence between the temperatures of the environment and 
of the superconducting leads. We believe that the low 
bias regime is experimentally most interesting because 
the junction in this case has low noise, low impedance 
and low thermal conductance, which opens up an oppor- 
tunity of its application as a sensitive detector. 

Apart from that, we study the noise of the heat cur- 
rent. It is also an important characteristic of the junction 
in view of its potential detector applications. 

On the technical side, we perform the usual second or- 
der perturbation theory in tunnel Hamiltonian connect- 
ing the two superconducting leads. In doing so we use 
the full counting statistics (FCS) approach to the heat 
transport which has been recently successfully applied to 
similar problem s 10 ' 11 . The method has an advantage of 
providing the information about the cooling power of the 
junction, its noise and, eventually, all cumulants of the 
heat current at one step. 

The paper is organized as follows: in Sec. II we intro- 
duce the model; in Sec. Ill we derive FCS of the junction; 
in Sec. IV we derive the general expression for the cool- 
ing power; in Sec V we consider the heat current noise; 
in Sec. VI we study average cooling power under various 
bias conditions; and, finally, in Sec. VII we summarize 
our results. Some details of the calculations are given in 
the appendices. 



II. MODEL 

We consider a Josephson tunnel junction with the nor- 
mal state resistance R connecting two superconducting 
leads with the energy gaps Ai, A2 and the temperatures 
TifT?. The junction is shunted by a capacitor C and by 
a resistor R$ kept at temperature T* , as shown in Fig. [T] 
The ohmic resistor provides the simplest model for the 
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FIG. 1: A Josephson junction between two superconductors 
with the gaps Ai, Ag and temperatures T\,Ti. The junction 
is shunted by a capacitor C, a resistor Rs and the whole 
circuit is biased by a current I x . The dynamics of this circuit 
is described by RCSJ model ©. 



electromagnetic environment of the junction. Our pri- 
mary goal is to find the heat current through the junction 
or, more precisely, the cooling power of the superconduc- 
tor 1 The latter is defined as the amount of energy 
extracted from the superconductor 1 per unit time. We 
will also be interested in the noise of the cooling power. 
The junction is described by the Hamiltonian 



H — H± + H 2 + Hj 



(1) 



where 



Hi = eifec; k c ak + Axc^kCi-k + Aicj_ fc cj 



Ik 



(2) 



is the Hamiltonian of the superconductor 1, 

H 2 = t2ka\ k a ak + A 2 a tk a^ k + A 2 a\_ k a\ k 

k L er=t,4. 



(3) 



is that of the superconductor 2, and 
Ht = V V \t kn e^ 2 al n c ak + t* kn e-^l 2 cl k a 



kn a— t,J, 



1,(4) 



is the tunnel Hamiltonian. Note that by means of a suit- 
able gauge transformation we made the superconduct- 
ing order parameters in the leads, Ai,A2, real. After 
this transformation the Josephson phase (p is moved into 
the phase of the tunneling amplitude. We will treat the 
phase ip as a quantum operator, which is indicated by 
the hat. We further assume that the quasiparticle distri- 
bution functions in the leads are Fermi functions, 



fj(E) = l/(l + e E ' k * T *), j = 1, 2. 



(5) 



Thus we ignore possible charge imbalance^ in the leads. 

Next, we will assume that the shunt resistance is much 
smaller than the resistance quantum R k = h/e 2 and the 
resistance of the junction 



Rs <C R k , R- 



(0) 



In this case, and at sufficiently high bias current and/or 
temperature, one can treat the Josephson phase as a 
classical variable. Its dynamics is then described by the 
RCSJ differential equation^ 



C 



Tup 
2e 



1 Kcp 
R~s^ 



Ic sin ip = I x + £s 



(7) 



where £s is the noise of the resistor Rs- At high temper- 
ature it is white Gaussian noise with the correlator 



Rs 



(8) 



The critical current of the junction, Ic , is determined by 
the integral^ 



^ dE _ A 1 A 2 [\-2f 1 {E)] 



eR J Al v /^_ Aj^Aj -£ 2 



(9) 



Here we have chosen Ai < A2, i.e. superconductor 1 
is supposed to be weaker than superconductor 2. 

Equation (|7|) describes the motion of an effective parti- 
cle with the coordinate (p in the tilted Josephson potential 



U{ip) = — Ej cos ip — hl x (p/2e, 



(10) 



where Ej — hlc/2e is the Josephson coupling energy. 
In order to ensure classical behavior of the phase in the 
wide range of parameters we also require 



Ej » E c , 



(11) 



where Ec — e 2 /2C is the charging energy. In what fol- 
lows we will focus on the regime of weak noise and require 
the condition 



k B T* < Ej 



(12) 



to be satisfied. This regime is the most interesting for a 
simple reason: if the environment temperature becomes 
comparable to or exceeds the Josephson energy, the su- 
perconducting branch of the current-voltage character- 
istics becomes strongly smeared and the heat transport 
properties of the Josephson junction become qualitatively 
similar to those of an NIS tunnel junction, which is well 
investigated. 

At low bias currents, I x < Ic, and at sufficiently low 
temperatures the classical approach based on Eq. Q is 
insufficient and the quantum nature of the phase becomes 
important. In this regime the phase is fluctuating around 
its equilibrium value 



ip = arcsin(J x /Ic), 



(13) 



and the charge and heat transport properties of the junc- 
tion may be expressed via the two quantum phase corre- 
lation functions^ 



V(E) 



,((v5(0)-^o) 2 )/2 



dt 
2nh 



e iEt/H / e i0(t)/2 e i0(O)/2 



(14) 
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These functions are expressed via the effective impedance 
of the electromagnetic environment of the junction, 
which in our model reads 



Z s (w) = [ -iwC + — + 2e f C cos<^o 
Ks —iWjj 



(15) 



The functions ([Til) may be explicitly evaluated and have 
the formi^ 



P(E) = 9(E rnax -\E\) J ^L t 

/ji 
2^ 6 



iEt/h+J(t) 



(16) 



iEt/H-J{t) 



(17) 



J(t) 



R s f dui oj (coth ^-(coswi — 1) — i sinwi) 



Rk 



+ R 2 C 2 {ufi-v*) 




where wj = ^2eIc/hC(l — I 2 /^) 1 / 4 is the frequency of 
small oscillations around the minimum of the Josephson 
potential. Here we have also introduced an effective high 
energy cutoff 



which equals to height of the barrier in the tilted Joseph- 
son potential pTJ|) . This cutoff comes from the fact that 
at higher energies the Gaussian approximation used to 
derive the correlation functions (|16I17[) is no longer valid. 
The high energy cutoff (|18[) is important only in the nar- 
row range of bias currents close to Ic ■ 



III. STATISTICS OF THE HEAT TRANSPORT 

In order to introduce the concept of the full count- 
ing statistics of the heat current we first define the en- 
ergy £\ extracted from the superconductor 1 during a 
time t. Since the heat transport through the junction is 
a stochastic processes, the energy £\ is a random value 
described by a certain probability distribution W{t,£\). 
Our goal in this section is to find the cumulant generating 
function corresponding to this distribution 



F{t, x) = In 



d£e- i£lX / h W(t,fi) 



(19) 



where \ is the counting field. The function (|19l) can be 
expressed in the form 



tr [ e - lH ix/fr e -iHt/h e iH 1 x/h e -H/k B T e iHt/h^ 



111 



tr [e 



-H/k B T] 



20) 



Provided the cumulant generating function is known, the 
cooling power is given by its derivative 



dt dt d\ 



(21) 



Similarly, zero frequency spectral density of the heat cur- 
rent noise is defined as follows 







d d 2 T 



dt 



dt dx 2 



(22) 



Both the cooling power and the noise fluctuate in time 
following the fluctuations of the phase ip. Hence, in order 
to obtain experimentally relevant parameters one should 
further average Eqs. (I21I22|) over these fluctuations. 

One can derive an explicit expression for the cumulant 
generating function (|2"0"|) by means of the second order 
perturbation theory in the tunnel Hamiltonian Ht- The 
details of this derivation are summarized in Appendix A. 
Here we directly cite the result. The function F(t,x) is 
the sum of two contributions: 



(23) 



where J- qp (t, x) and J~,j{t, x) originate from, respectively, 
quasiparticle and Cooper pair tunneling through the 
junction. Defining the two functions characterizing the 
superconducting leads (j = 1,2) 



WP(t) - 



dc j 



where Ej = J e 2 + A 2 , and, assuming that the Joseph- 
son phase is classical, we obtain the quasiparticle contri- 
bution in the form 



' ?*>f dt"^~^ 



■"• ' -h< -l< ./, - ./ , - 2 
x {W? p (t' -t" + x )Wf {t'-t") 

+ w? p (t" -t' + x )w 2 qp (t"-t')}. 



(25) 



The Cooper pair, or Josephson contribution has a similar 
structure, 



Tj = 



1 



dt' 



t' 



dt" cos 



¥>(*') + ¥>(*" 



■Khe 2 Rj Q ' 7-oo " 2 
x{W 1 J (t'-t" + x )W 2 J (t'~t") 
+ W 1 J (t"-t' + x,)W 2 J (t"-t')}, (26) 

but it depends on the sum of the two phases, (p(t') + 
<p(t"), instead of their difference. Besides that, it contains 
another pair of the functions 



wHt) 



de i 



'2Ei 



x [e^ME,) - e- iE ^ n (l - ME,))] . (27) 



The expressions (|25I26[) are valid for arbitrary time de- 
pendence of the Josephson phase (p{t). 

If the phase is quantum, the cumulant generating func- 
tions have to be modified and take the form 



qP ~ 2^R 



dt' / dt'' 
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x I ^ e -i0(t')/2 e i#(i")/2^ + ^ e <0(t')/2 e -i0(t")/2^ 

xiy 1 9P (i'-t" + x)w r l p (t , -0 

+ n e -i<p(t")/2 e iip(t')/2^ + ^ e #(t")/2 e -i#(t')/2\^ 

X (*" - t' + X )Wf (t" - o), (28) 



1 



2whe 2 R 







dt' / di'' 



< ^e^( t ')/ 2 e ^( t ")/2^ _|_ ^ e -i#(t')/2 e -i#(t")/2^ 

xTy 1 J (t'-t" + X )^ 2 7 (<'-0 

+ ^/ e ^(*")/2 e »V(*')/2^ + ^ e -#(t")/2 e -iv(*')/2\^ 

xW/(t"-t' + X )W/(i" -*')}• (29) 



The angular brackets here stand for both quantum me- 
chanical and statistical averaging. At bias currents below 
Ic the fluctuations of the phase become Gaussian, and 
the average phase correlators, appearing in Eqs. (|28I29I) . 
may be expressed via the functions (|16ll7p . At I x > Ic 
the classical dynamics of the phase sets in, and under the 
conditions (|6|11|12|) one can replace the cumulant gener- 
ating functions (|28I29[) by the classical ones (I25I26P . 



IV. COOLING POWER 

We are now in position to evaluate the cooling power 
of superconductor 1. In the regime of classical phase 
fluctuations Eqs. (|21I25I26[) lead to the following result 



... ■ ■ 1 ■ ./ ■ (30) 
where the quasiparticle cooling power has the form 

p &Ht) = — b f dt '[w? p (t - t'wnt - f ) 

7T6 it J _oq 

+ Wf{t' - t)Wf(t! - t)] cos ~ , (31) 
and the Josephson contribution reads 



ire 2 R 



dt'[W/(t-t')Wi{t-t') 



+ Wftf - t)W 2 J (t' - t)] cos y(t)+ / (f,> . (32) 



It is interesting to compare the results (|30I31I32[) with 
the well known expressions for the charge curren t 15 ! 16 



/(<) = I gp (t) + Ij{t), 



(33) 



nheR 



dt'{W? p (t-t')W2 P (t-t') 



Wf(t' - t)Wl P (t' - t)} sin Vit) ^ , (34) 



Ij(t) = 



nheR 



df {w/ (t - t')wi \t - 1') 



Wl J (t'-t)W 2 J (t'-t)} 



(35) 



One observes two main differences between the cooling 
power and the charge current: (i) the charge current 
contains sines of the combinations ip(t) ± (p(t') while the 
cooling power - their cosines, and (ii) the charge current 
contains the functions W^ p , W( while the cooling power 
— their time derivatives W^ p , W( . 

Combining Eqs. ()31I32|) with (|34I35[) one can easily 
derive the identity reflecting the energy conservation in 
the junction 

pa> + P» = !g (w „) + ^, (36) 

where P^ is the cooling power of superconductor 2, 
which differs from PW (|30fl32|) by interchanging the in- 
dexes 1 and 2, and A is a certain combination of the in- 
tegrals, which vanishes if the phase varies in time slowly. 
The last term on the right hand side of Eq. (f36| drops 
out after the averaging because it is a full time deriva- 
tive. Obviously, the product hipl qp /2e is the work done 
by the external current source. The term Tuplj/2e in 
the adiabatic limit reduces to the time derivative of the 
Josephson energy —d(Ej cos ip) / dt and vanishes upon the 
averaging. Hence the identity (|36| shows that the work 
of the current source is, on average, split between the two 
superconducting leads. Moreover, if the lead 1 is cooled 
and PW > 0, then according to Eq. (|3"o]l lead 2 is in- 
evitably strongly heated because it has to absorb both 
the power pW and the work of the current source. This 
effect is well known for NIS tunnel junctions 1 . 

Let us now consider a voltage biased Josephson junc- 
tion and put ip(t) = 2eVt/h. The cooling power takes 
the form 

P«(i) = Pg\V)+P^(V)co S [2eVt/h] 

+ P^(V)sm[2eVt/h}. (37) 

The quasiparticle contribution to it, Pq P \v), becomes 
time independent in this case and is given by the integral 5 

e"2p, 

x(E- eV) [h(E - eV) - f 2 (E)] . (38) 



p^Hv) 



, dE N X {E - eV)N 2 {E) 



Here we defined the quasiparticle densities of states in 
the leads 



Nj(E) = \E\9(\E\ — Aj)/ J E 2 — Aj, j = l,2. (39) 
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FIG. 2: The quasiparticle cooling power P qP (V) (solid line) 
and the absolute value of the amplitude |P CO s(V)| (dotted 
line). The parameters are chosen as follows: Ai = 100 ^teV, 
A 2 = 200 iieV, R = 1 kfi, Ti = T 2 = 230 mK. 




FIG. 3: The Josephson component of the cooling power 
P s in(V) for the same parameters as in Fig. [5] 



The amplitude of the " anomalous" oscillating term reads 



1 

AiA 2 
E 



dE Ni(E - eV)N 2 (E) 
[h(E - eV) - f 2 (E)], (40) 



and the Josephson amplitude has the form 

Mr- V 



2ireR 
1 - fi{Ei) - h(E 2 ) 



a a A i A 2 
de\de 2 — — — x 



+ 



E 2 

h(E 1 )-h{E 2 ) 



{E 1 +E 2 ) 2 -e 2 V 2 (E 1 - E 2 f - e 2 V 2 



.(41) 



The amplitudes (V) and P^os (V) are exponentially 
suppressed in the limit of low bias voltage and low tem- 
peratures eV, Ti, T 2 <C A2 — Ai- In this case we find 



P^(V) = 



5/2 



2 R VA 



2 _ A 2 
2 ^1 



x cosh ■ 



eV 



y/k^F 2 e- A */ k ° T > 



k B Ti 

P^l(V) = -(A 1 /A 2 )P gp (V). 



(42) 



In contrast to that, the Josephson amplitude P s \ n ' (V) 
is not exponentially suppressed. However, it is a linear 
function of V . 



pM(V) = a(A 2 /A 1 )I c V, 



(43) 



and tends to zero at V = for any values of the tempera- 
tures Ti, T 2 . Here we have defined a numerical pre- factor 



1 



K' (VT^ 2 ) 



vi - k 2 A"(vr 



(44) 



where if(ir) is the complete elliptic integral of the first 
kind. 

In the limit of zero bias, V — > 0, the cooling power (|37|) 
reduces to the former- 



pa) = _ 
e 2 i? 



E(E 2 - AiA 2 cos<^o) 
V^-A 2 /^!! 



x [h{E)-f 2 {E% 



(45) 



and only the contributions Pqj) and Pioi survive. 

If the voltage drop across the junction is not 
fixed, but the phase changes adiabatically, 99 <C 
min{Ai, A2, Ti, T2}, one can approximate the cooling 
power by replacing V by h(p/2e and 2eVt/h by if in Eq. 
(|37p . After that one arrives at the approximate expres- 
sion 



P«(t) 



p(i) «£ 
9P I 2e 



2 

A 2 



, P (i) f'M 
2e 



COS (y5 



(46) 



The last Josephson term of this formula is the full time 
derivative and it averages out to zero for any realization 
of the phase fluctuations. Due to the same reason the 
Josephson component results in the contribution cx uj 2 to 
heat current noise, which vanishes in zero frequency limit. 
In order to detect it one has to perform high frequency 
noise measurements. 

The components Pqp and Pioi of the cooling power 



are plotted in Fig. [5] and the Josephson component P t 



(i) 



is shown in Fig. [3l We observe that Pq P ' and P CO s are 
even functions of the bias voltage, while the component 
Pain ^ s an 00 -d function of it. 



V. HEAT CURRENT NOISE 

Low frequency heat current noise, or the noise of the 
cooling power, S P = J dt' (SP ^(t)SP^(t')), can be de- 
rived from the Eqs. (|22I25I26[) and has the form 



qqp 
Op 



S p . 



(47) 
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Here 



ClP 
ij p 



Ri 



2n 2 R 

+W? p {t' -t)Wf{t' -*)} 



dt'{W? p {t-t')Wf(t-t') 

<p(t) - ip{t>) 



(48) 



is the noise associated with quasiparticle tunneling and 



S\ 



2tt 2 R 

W?{t' -t)W£{t' -t)} cos 



dt' {Wf (t - t')Wi (t - t') 

<p(t) + ¥>(f ) 



(49) 



is the Josephson contribution to the noise. 

Considering again a voltage biased junction with 
(p(t) = 2eVt/h, we find±I 

Sf = -^JdEN 1 {E~eV)N 2 (E)(E-eV) 2 

x [ME - eV)(l - h(E)) 

+ (l-f 1 (E-eV))f 2 (E)]. (50) 

The Josephson noise in this case takes the from 

Sp = S cos {V) cos[2eVt/h] + S ain (V) sm[2eVt/fi}, (51) 

where 



S cos (V) 



1 



dE Nx(E - eV)N 2 {E) 



e 2 R 
A^jE-eV) 
' E 



and 



S ai n(V) 



+ (l-f 1 (E-eV))f 2 (E)], 

A 1 A 2 E 1 



[fi(E - eV)(l - h(E)) 
(52) 



V 



de\de 2 - 



2ireR J E 2 
MEi)ME2) + [1 - MEiW - ME 2 )} 



{E!+E 2 ) 2 



2y2 



MEi)[l - ME2)} + [1 - MEi)]ME2) 



{Ei - E 2 ) 



2y2 



(53) 



Scos{V) and 5 s i n (V) are the amplitudes of oscillating 
components of the noise and unlike S qp (V), which is al- 
ways positive, may change their sign at certain values of 
the bias voltage. For example, the amplitude S cos (V) is 
negative at e\V\ < Ai + A 2 and positive otherwise. 

In a typical experimental setup the time averaged value 
of the noise is measured. In this case the oscillating com- 
ponents vanish and the noise spectral density is given by 
the quasiparticle noise Sff (|50)l . 

At zero bias voltage, V = 0, the phase takes its equi- 
librium value ifo (|13[) and we get 

S P = ^ J dEN 1 (E)N 2 (E)(E 2 -A 1 A 2 cos(p ) 

x [ME)(l - ME)) + (1 - ME))ME)] • (54) 
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FIG. 4: The components of the heat current noise S q p(V) 
(solid line) and the |S , ™ S (V)| (dotted line). The parameters 
of the junction are the same as in Fig. ([2} • 



One can verify that in equilibrium, i.e. at T\ = T 2 = T, 
zero bias noise satisfies the fluctuation dissipation theo- 
rem in the form 



S P = 2k B T 2 



dpW 



9Ti 



(55) 



T 1= T 2 =T 



where is given by Eq. (|45|) . 

The noise for typical values of the junction parameters 
is plotted in Fig. 0] It exhibits characteristic peaks at 
bias voltages eV = ±(A 2 — Ai) and jumps at eV^ = 
±(A 2 + Ai), which are also seen in the cooling power 
and the current. 



VI. EFFECT OF PHASE DYNAMICS 

So far we have discussed a voltage biased Josephson 
junction and a junction in superconducting state with a 
fixed phase. In practice these regimes are not easy to 
achieve and, in general, one has to consider a junction 
in combination with its electromagnetic environment as 
depicted in Fig. [TJ The effect of the environment, and 
of the phase dynamics associated with it, on the cooling 
power of the junction will be the subject of this section. 

We begin with the regime of high bias current and 
assume that the conditions (|6I11I12|) are satisfied. The 
phase in this case may be treated classically and obeys 
Eq. ([7]). In order to simplify our analysis further, we 
will consider an overdamped Josepshon junction assum- 
ing that 1/RsC ^ \/2eIc /fiC '. This condition allows 
one to put C = in the equation of motion ([7]), which 
reduces it to the exactly solvable^ resistively shunted 
Josephson junction model (RS J- model). Some details on 
the RSJ model needed for our analysis are briefly summa- 
rized in Appendix B. Formally, in order to evaluate the 
average value of the cooling power one needs to resolve 
Eq. (JZl) for ip(t), substitute the result in Eqs. pi|32p and 
evaluate the path integral over the noise of the shunt re- 
sistor £,s(t)- However, at high values of bias current one 
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can avoid doing that and instead use a simple approxi- 
mate approach. In order to justify it, we note that the 
integrals in Eqs. (|31I32|) converge if t - t' > h/2eV{t), 
where V(t) — %ip(t)/2e is the instantaneous value of the 
voltage. If the condition 



i([M* + */2eV(t))-M*)] a 



eR s k B T* 
2hV(t) 



< 1(56) 



is satisfied, one can ignore the noise in Eq. <[7|) while 
solving it on the short time intervals of the duration ~ 
h/2eV(t). Then from Eq. © one finds 

V(t) = R s {I x -I C sin<p(t)), (57) 

and the voltage is restricted to the interval 

Rs(I x - Ic) < V{t) < R S (I X + I c ). (58) 

The condition (|56[) then translates into a simpler one 

h ~ Ic » \k B T* . (59) 

Thus, provided this condition is satisfied, one can make 
the replacement (|57j) and perform the noise averaging in 



Eqs. (|30H35|) with the aid of the following approximate 
formula 



,i[tp{t)-<p{t')]/2 



dip e 



ieR s (Ix-Ic sin ip)(t-t')/h 



a(ip). 



(60) 



Here a(<p)12- is the distribution function of the phase (|B1[) . 
As a result, the average current (|33I34I35[) and the cool- 
ing power (|30I31I32[) acquire the form 



(I) 
( p(U) 



/7T 
dipa(ip) [l qp (I x Rs - IcRs sin p) 
-7T 

+ I c sm(p], (61) 
/ dipa(ip) [P qp (I x R s - IcRs sin <p) 

J — 7T 

+ -Pcos (4 Rs - IcRs sin ip) cos tp], (62) 



where 



I gP (V) 



dE 



N X {E - eV)N 2 (E) 
eR 



[f x {E-eV)-h{Ej\ 



is the usual quasiparticle current. 

The current (|6"Tj) and the cooling power (|62|) are plot- 
ted in Fig. [5] versus the average voltage drop across the 
junction (V) (|B2|) . As expected, strong noise of the shunt 
resistor suppresses the superconducting branch of the I- 
V curve at high temperatures. At the same time, it does 
not drastically change the cooling power as compared to 
the voltage biased case shown in Fig. O The only effect 
of the noise at high bias (15^1) is the smearing of the peaks 
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FIG. 5: The time averaged current (/) (|61[l (top panel) and 
cooling power (P) (|62p (lower panel) at different temper- 
atures. The parameters are the same as in Fig. [2l i.e. 
Ai(0) = 100 ^teV, A 2 (0) = 200 fj>eV, R=lkQ. The tempera- 
ture dependence of the energy gaps Ai(T), A 2 (T) is modeled 
by Bardeen-Cooper-Schrieffer theory. Besides that we have 
put Ti = T 2 = T* and chosen R s = 100 Q.. The horizontal 
axis in both plots shows the average voltage drop across the 
junction (V) given by Eq. (|B2[) . 



located at eV = A2 — Ai. This smearing is caused by 
the fluctuations of the voltage in the interval ([55]) . 

Next, we consider a simple example of an overdamped 
junction shunted by a noiseless resistor and put £<j = 
and C — in the dynamical equation ([?])• This regime 
is formally achieved in the limit 1/RsC >• \j2elc /hC 
k B T*/Ej — > and e 2 Rs/h -> 0. The model becomes 
exactly solvable in this case because the condition (|5"9")l 
is always satisfied at I x > Ic, and at I x < Ic the phase 
is just pinned to the equilibrium value tpo (fT3|) . Hence 
at I x < Ic the cooling power is given by Eq. (|45|) . At 
Ix > Ic the Eq. ([7]) can be solved exactly and one finds 



V(t) 



I 2 c)Rs 



2e 



Ix - Ic sin u t 



(63) 



Here u>o — 2eRs\/I x — I c /h is the Josephson frequency 
proportional to the average voltage (V) — Rs\/T!f 
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FIG. 6: The average cooling power {P^ 1 ') (|45l64p of a re- 
sistively shunted Josephson junction at low bias and in the 
absence of the noise. Top panel: (P ) as a function of the 
bias current at T\ = T 2 , Ti > T 2 and Ti < T 2 . Lower 
panel: (-P' 1 ') as a function of the time averaged voltage 
(V) = Rs y/ Ix — Iq\ for comparison we have also shown 
the quasiparticle cooling power P qp (|42p of a voltage biased 
Josephson junction. The parameters are chosen as follows: 
Ai = 100 fieV, A 2 = 200 ^eV, R = 1 kfi, Rs = 0.1 kQ. For 
these parameters we find Ic = 212 nA and IcRs = 21.2 //V. 



The average cooling power is found by the time averaging 
of the Eq. (l4"6l and reads 

(P(D) = JL / dt{P$(y)+co SV >P£l(V)}. (64) 

One can verify that in the limit fcgT* /Ej — > the 
Eq. (|64| matches Eq. (l62t because in this limit the 
phase distribution function (|B1[) reduces to the form 
a ((f) = \Jll~ I c /2ir(I x - I c sirup) for I x > I c and 
a ((p) =S(tp — arcsin( I x /Ic)) for 4 < ^c- 

The time dependence of the cos ip at I x > Ic reads 



costp = y J2 - /2 cosuj t/(I x — Ic sinwo*)- (65) 

With the aid of this result one can show that the term 
oc P cos in Eq. ([(14]) identically vanishes. 



In the most interesting low bias and low temperature 
regime, eV < k B Ti,T 2 <C Ai,A2, we can approximate 
the quasiparticle cooling power by a simple formula which 
follows from the Eq. (g2J), 




(66) 



2 e 



k B T 2 



Afterwards, the time integral in Eq. (|64p is easily evalu- 
ated, and at I x > Ic one arrives at the result 



(P 



= 5P((V)) 



A 5/2 
27T A 2 



^k B Ti e fc B^i _ v / k B T 2 e k B^ 



e 2 R 



Here the contribution 



a/at 



A? 



{.67) 



FAfe- A ^^ T i|(y)| Vt'cR's + ( v ) 2 , 



2 (fe B Tx)3/VAi - A? 



i? 



-(68) 



arises from voltage oscillations (155)) . 

It is clear from Eq. (|6"51) that at low voltage, (V) % 
IcRs the cooling power grows linearly with the voltage, 
SP oc |(V)|. This unusual behavior may be understood if 
one recalls that phase slips are responsible for the trans- 
port in this regime. Indeed, at (V) < IcRs the time 
dependence of the voltage (|6"5)) reduces to a series of well 
separated pulses corresponding to a sequence of phase 
slips. Hence one can express the average cooling power as 
a product of the phase slip rate, which is robustly related 
to the average voltage as r ps = e|(V)|/7wr by Josephson 
relation, and the energy carried away from superconduc- 
tor 1 by a single phase slip E ps , 



SP((V)) = E ps T ps . 



Comparing equations (|68|) and (|69|) we find 



E„ 



P$(V(t))dt 



-7r/W0 



r 3 / 2 ftA^V^/^ v/7p| + (V)> 



V2 {k B TxYt*y/E 

7T 3 / 2 R S 



A? 

A 2 /k B T 1 



eR 



K\I X 



V2 R (fc s Ti)3/VA3 - A 2 



(69) 



(70) 



(71) 



(72) 



The energy E ps becomes voltage independent for (V) < 
IcRs because the phase slips become well separated. 
Hence Eq. (|f]^|) always leads to the linear dependence 
of the cooling power on voltage in this limit. This con- 
clusion holds for more general models of Josephson dy- 
namics. For example, weak noise, k B T* <C Ej, or finite 
junction capacitance do not change it. 
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The dependence of the cooling power of a noiseless 
Josephson junction on the bias current I x and average 
voltage (V) is shown in Fig. El 

In the remaining part of this section we will gener- 
alize the model introduced above and include the weak 
noise with fc^T* < Ej into the analysis. Two types 
of phase trajectories may be distinguished in this case: 
phase slips and Gaussian fluctuations around a minimum 
of the Josephson potential U((p) (|T0|) . Accordingly, the 
cooling power may be expressed as a sum of two terms 



(P (1) ) 



P 



ps 



P 



gs- 



(73) 



where P ps is the contribution of the phase slips and P gs 
is the contribution of the Gaussian fluctuations. 

Let us first consider the phase slip contribution. As we 
have discussed above, it is given by the product 



P = ^E 



(74) 



In the presence of noise the average voltage is given by 
Eq. (IB2|) while the phase slip energy at \I X \ > Ic is still 
defined by Eqs. (JTT|) or ([72]). At \I X \ < I c every phase 
slip is accompanied by a voltage pulse of the form 



V ps (t) 



Rs{I 2 c 



/ccosh(2ei?sV^T 



lit in 



(75) 



which may be derived for Eq. ([7]) with C = and £s = 
0. The corresponding phase slip trajectory, ip ps (t) — 

J* dt'2eV ps (t') /%, connects the maximum, — it — p , and 
the minimum, ip a , of the Josephson potential. The phase 
slip cooling energy then takes the form 



dtP$(V ps (t)) 

) 

7T hR s A 5 2 /2 e - A -/ kBTl 



2 R{k B T x f/ VA1 - A? 
21 



I 2 



arctan 



Ic 



In -I,. 



(76) 



which exactly matches the expression (1721) at I x = Ic- 

Let us now derive the contribution of Gausian fluc- 
tuations P gs . Provided the subgap quasiparticle resis- 
tance of the junction exceeds the quantum resistance, 
ft e A 2 /k B T 1 s^, j^ k ^ w hich is often the case in experiments 
with small area junctions, we may evaluate P gs apply- 
ing the theory of dynamical Coulomb blockade 14 . In this 
way we arrive at the result 



Pgs = ^RJ dE i dE M E i) N 2( E i) E i 
x {fi(E 1 )[l-f 2 (E 2 )]V(E 1 -E 2 ) 
-[I- h(E 1 )]f 2 (E 2 )V(E 2 -E 1 )} 



,-(V>/2 



COS Ifio 



5 2 R 



dE 1 dE 2 N 1 {E 1 )N 2 {E 2 ) 



^^{A^OIl - f 2 (E 2 )]V(E l - E 2 ) 
[1- h(E 1 )]f 2 (E 2 )V(E 2 -E l )}. 



(77) 



Here the we have introduced the average square of the 
phase fluctuations 



e 2 Rs 



duj- 



uj coth tjt^t 

2k B l 



uj 2 + R 2 C 2 {lu 2 - ui 2 )' 



(78) 



and the functions V(E), V(E) are defined in Eqs. (|16I17I) . 
They obey the detailed balance principle 

V{-E) = e E ' kBT *V{E), V{-E) = e E / kBT "P(E),{79) 

which ensures that P gs = at T\ = T 2 = T* , as required 
by thermodynamics. We also note that in the limit R$ — > 
one finds V(E) = V(E) = 5(E), (Sip 2 ) = and the 
expression (1771) reduces to Eq. (|4"5"|) as expected. 

In the limit Rs <C Rk one can expand the function 
V(E) (|16p in powers of J(t), which leads to the approxi- 
mate expression 



V(E) 



e 2 Rs 



0(E n 



\E\ 



^ E 2 + R 2 S C 2 (E 2 — h 2 uj 2 ) 
E (Sip 2 ) 



1 _ e -E/k B T* 



S(E).(80) 



The function V(E) in this limit differs from V(E) only by 
the sign in front of the first term. Within this approxima- 
tion and for equal temperatures of the superconducting 
leads, Ti = T 2 = T, the Gaussian contribution to the 
cooling power (jTT[) takes the form 



Pa. 



2e 2 R s 



X P r , 









t 


dE 














i 




E 2 


f R 2 C 2 (E 2 - 



2 COS ipo 



E 



E 



oE/k B T* _ l e E/k B T _ l 



(81) 



The average cooling power affected by the noise 
([73I74I81P is plotted in Fig. [7J Comparing it with the 
cooling power of the noiseless junction, Fig. [BJ we ob- 
serve two major differences. The first one is an expected 
smearing of the sharp features of the (P)(I X ) dependence 
by the noise. The second difference becomes visible at 
\Ix\ < Ic an d T\ = T 2 = T, where the cooling power 
vanishes in the absence of the noise but remains finite 
if the noise is present. The noise induced cooling power 
is positive if T* > T and negative in the opposite case 
T* < T. Thus, even at zero bias current one can cool 
the superconductor 1 by heating the shunt resistor and 
heat it by cooling the resistor down, see the inset in the 
top panel of Fig. [7] This interesting effect has been first 
predicted for an NIS tunnel junction 19 . 
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FIG. 7: Time averaged cooling power (P) (|73l74l8ip at Ti = 
T2 — 230 mK and at various values of the noise temperature 
T* . Other parameters are the same as in Fig. [6] Top panel: 
Cooling power as a function of the bias current, inset: zero 
bias cooling power versus the environment temperature T*. 
Bottom panel: Cooling power as a function of the function of 
the average voltage for three different values of T* . 



Finally we briefly note that although phase dynamics 
significantly modifies the average cooling power at low 
bias, it only weakly affects the heat current noise. Essen- 
tially, one can always use the equations (I50I54[) to make 
estimates of the noise. The reasons for that are simple. 
First, the energy carried by one quasiparticle and the 
rate of quasiparticle tunneling in most cases exceed, re- 
spectively, the energy carried by one phase slip and the 
rate of phase slip tunneling. Second, the noise remains 
finite even at zero bias and at T\ = T2 = T* . Thus 
its averaging over phase motion may only lead to a cer- 
tain smearing of the bias dependence of the noise without 
qualitatively changing it. 



VII. SUMMARY 

In summary, we have analyzed the heat transport 
through a Josephson tunnel junction. We have devel- 



oped the full counting statistics approach to this prob- 
lem and derived general expressions for the cumulant 
generating function of the heat extracted from one of 
the superconducting leads (|23H26[) , for the corresponding 
cooling power (|30ti32[) and for the heat current noise (|47l - 
I4U1) . These general expressions are valid for an arbitrary 
time dependence of the Josephson phase. In analogy 
with the charge current, the cooling power is the sum 
of the quasiparticle and Josephson contributions. The 
latter contribution oscillates in time and averages to zero 
if the junction is biased above the critical current. We 
have also generalized all these results to the regime where 
the quantum nature of the Josephson phase becomes im- 
portant. Combining these findings with the conventional 
theory of the Josephson phase dynamics, we have derived 
the average cooling power of the junction, which can be 
measured in the experiment. We find, in particular, that 
at low bias the average cooling power is proportional to 
the voltage drop on the junction. Finally, we have also 
shown that zero bias cooling power of the junction may 
be tuned by the heating or cooling the resistor connected 
in parallel to it. 

We acknowledge support from . . . 



Appendix A: Explicit expression for the generating 
function 



In this appendix we derive an explicit expression for 
the cumulant generating function Fit, \) (pDl by means 
of the perturbation theory in tunnel Hamiltonian. The 
first step of our derivation is to transform F{t, \) to the 
form 



J-(i, X )=ln 



tr y e -iH x t/h e ~H/k B T e iHt/h^ 



(Al) 



/tr [e-v/kBT} 

where we have introduced the transformed Hamiltonian 

(A2) 



H x = e - iHlx/n He iHlx/h . 



Since the Hamiltonians H\ and H2 commute with H±, 
only the tunnel Hamiltonian changes under this trans- 
formation. Thus we get 

H x = H 1 +H 2 + H T (x), (A3) 

where 

H T ( X ) = e - lHlX ' h H T e lHlx/h 



£T kn 



(A4) 



The operators cvfc(x), c^ fc (x) in this expression are de- 
fined as follows 



Qxfc(x) 

4« 



_ p -iHi_x/n 



iH lX /h 



e -iH lX /n c t e iH lX /h_ 
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Next we perform Bogolubov transformation of the op- 
erators: 

Cffc = Mlfe7tfc ~ v ^l\^ k , Clk = vikll^ k + u ^lik, 
a^k = u 2k a^ k - V2kOi\_ k , a^k = v 2k a\_ k + u 2k a^ k {.P^) 

Here 7, 7^ are the quasiparticle annihilation and creation 
operators in lead 1 and a, are similar operators acting 
in lead 2. Besides that we have introduced the coherence 



factors (here Ej k = J e 



-jk 



Ujk 



V2 



e jk 



E 



jk 



Vjk 



1 
71 



1 - 



e jk 



E 



(A6) 



The averaging of the products of the quasiparticle oper- 
ators results in the distribution functions in the leads 

Uklak) = fi(E lk ), (aLa an ) = f 2 (E 2k ). (A7) 

After the transformation (|A5|) the Hamiltonian Hi ac- 
quires the diagonal form 

Hi = J>i fc - Ei k ) + Elk Hk-rtk + llklik] ■ (A8) 



The next step is two switch to the Heisenberg represen- 
tation. In order to do that we split the total Hamiltonian 
into two parts: the exactly solvable part Ho = Hi + H 2 , 
which includes only the Hamiltonians of the leads, and 
the perturbation V x = Ht(x)- Next we introduce the 
time ordered exponent S x (t), 



-W x t/h _ e — iH ot/hg ^ 



S x (t) — T exp 



dt'V x (t') 



where we defined the time dependent tunnel Hamiltonian 
V x (t) = e lHot/n H T ( X )e- lHot/n 



cr kn 
t *kn C 



- im/2 'Ut,X)a r rn(l)_- 



(A9) 



The time dependent creation and annihilation operators 
in this expression read 



ctk(t,x) 
Cik(t,x) 



u 2n e- iE ^ h a tk - v 2n e iE ^ n ^ 



J,— 71' 



-iE 2 „t/h 



Oi\.k + v 2n e 



u 2n e 



ui k e %Elkx/n e~ lElkt/h ^ k 

u lk e iE ^/ h e- iE ^ h %k 
+ vi k e- lE ^ h e lE ^ t ' h 1 \_ k . 



(A10) 



After these transformations the cumulant generating 
function (|A1|) takes the form 



E(t, X )=ln 



tr 



S x (t)e- H / k * T Sl(t)\ /tr \e- H ' k * T 
Performing a formal expansion of the operator S x it) in 
this expression up to the second order in V x we find 



dt< / dt"[- (V x (t')V x (t")) - (V (t")V (t')) 



(All) 



(V (t")V x (t')) + {V (t')V x (t"))]. 



It is now straightforward to evaluate this function. One 
should use the definition of the interaction potential V x (t) 
(|A9|) , combine it with the explicit expressions for the elec- 
tron creation and annihilation operators (|A10[) and apply 
Wick's theorem together with Eqs. (|A7|) to evaluate the 
average values of the products of these operators. Omit- 
ting the terms which do not depend on the counting filed 
X, one arrives at the results (|23|28|29|) . 



Appendix B: Dynamics of Josephson phase 

Here we briefly summarize well known results^ on the 
phase dynamics in a resistively shunted Josephson junc- 
tion, which is described by Eq. ([7]) with C = 0. The 
probability distribution function of the phase, <j(t,ip), 
satisfies the Fokker-Planck equation 



da 2eR s d 



dt h dip 

Ae 2 R s k B T* d 2 a 



[(4 ~ I c sin ip)a] 
= 0. 



dip 2 



The stationary solution of this equation, which is periodic 
in phase and to which any solution approaches in the long 
time limit, reads 



a(ip) 



h(V) 



4ireR s k B T* 



dip e 



lhI x (<p-<p')/2ek B T* 



x e 



(cos ip— cos ip')Ej /ksT* 



(Bl) 



This solution is normalized as follows: dipa((p) = 1. 
In Eq. (jBll) we have introduced the average voltage drop 
on the junction (V), which reads 



(V) = 



ksT_eRs. ginh kM^_ 

n 2ek b i 



;; /2 ^cosh(^)/ (^ 



(B2) 
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